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Abstract. We prove a priori estimates for the three-dimensional compressible Euler 
equations with moving physical vacuum boundary, with an equation of state given by 
p{p) = Cjp'^ for 7 > 1. The vacuum condition necessitates the vanishing of the pressure, 
and hence density, on the dynamic boundary, which creates a degenerate and charac- 
teristic hyperbolic free-boundary system to which standard methods of symmetrizable 
hyperbolic equations cannot be applied. 
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1. Introduction 

1.1. The compressible Euler equations in Eulerian variables. For < t < T, the evo- 
lution of a three-dimensional compressible gas moving inside of a dynamic vacuum boundary 
is modeled by the one-phase compressible Euler equations: 

p[ut + u ■ Du] 
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The open, bounded subset ri(t) C M.^^ denotes the changing volume occupied by the gas, 
r(<) := dn{t) denotes the moving vacuum boundary, V(r(t)) denotes normal velocity of 
r(i), and n denotes the exterior unit normal vector to r(<). The vector-field u = (ui, U2, u^) 
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denotes the Eulerian velocity field, p denotes the pressure function, and p denotes the density 
of the gas. The equation of state p{p) is given by 

p{x,t) = C^p{x,ty for 7>1, (1.2) 

where is the adiabatic constant which we set to unity, and 

p > in D.{t) and p = {) on T{t) . 

Equation (|l.lap is the conservation of momentum; (jl.lbp is the conservation of mass; the 
boundary condition (|l.lc[) states that pressure (and hence density) vanish along the vacuum 
boundary; (jl.ldp states that the vacuum boundary is moving with the normal component 
of the fluid velocity, and ()l.le|) - (ll.lf|) arc the initial conditions for the density, velocity, and 
domain. Using the equation of state ()1.2|) . ()l.lap is written as 

p[ut + u- Du]+ Dp'' m. Vl{t). (fTTal ) 



1.2. Physical vacuum. With the sound speed given by c := ^Jdp/dp and N denoting the 
outward unit normal to F, satisfaction of the condition 

§ < on r (1.3) 

defines a physical vacuum boundary (see [TU], [T^], [T3], [H], [H], [10])) where cq = c\t=o- 
The physical vacuum condition (|1.3p is equivalent to the requirement that 

M_!<Oonr. (1.4) 

Since po > in fJ, \1A\ implies that for some positive constant C and a; G near the 
vacuum boundary F, 

Po~^(2;) > Cdist(a;,F). (1.5) 

Because of condition (jl.Sp , the compressible Euler system (jl.ip is a degenerate and charac- 
teristic hyperbolic system to which standard methods of symmetric hyperbolic conservation 
laws cannot be applied. 

We note that by choosing a lower-bound with a faster rate of degeneracy such as, for 
example, dist(2;, F(i))'' for 6 = 2, 3, the analysis becomes significantly easier; for instance, 

Dp'^^^ (x t) 

if 6 = 2, then . " _^ is bounded for all x ^ fl. This bound makes it possible to 

VPo (^.*) 

easily control error terms in energy estimates, and in effect removes the singular behavior 
associated with the physical vacuum condition (jl.Sp . 

1.3. Fixing the domain and the Lagrangian variables on H.. We transform the system 
(jl.ip into Lagrangian variables. We let rj{x, t) denote the "position" of the gas particle x at 
time t. Thus, 

dfi] = uo ri for < > and ri{x, 0) = x 
where o denotes composition so that [u o 77] (x, t) u{ri{x, t),t) . We set 
V = uo rj (Lagrangian velocity), 
f — p ° rj (Lagrangian density), 
A = [Dr]]~^ (inverse of deformation tensor), 
J — det Dr] (Jacobian determinant), 
a = J A (tranpose of cofactor matrix). 
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Using Einstein's summation convention defined in Section [2. SI below. and using the notation 
F,k to denote ■§§-, the ki/i-partial derivative of F for k — 1,2, 3, the Lagrangian version of 
equations (ll.la|) - (ll.lb|) can be written on the fixed reference domain as 

fvl + Alr,k^O mnx{0,T], (1.6a) 

ft + fAlv\j^O mnx{0,T], {1.6b) 

/ = in X (0, T] , (1.6c) 

{f,v,ij) ^ {po,uo,e) mnx{t^O}, (1.6d) 

where e{x) = x denotes the identity map on fi. 

Since Jt — JAlv^,j and since J(0) = 1 (since we have taken r]{x, 0) — x), it follows that 

f^PoJ-\ (1.7) 

so that the initial density function po can be viewed as a parameter in the Euler equations. 
Let r := 9f2 denote the initial vacuum boundary; using, that A'^ — J^^ a\, we write the 
compressible Euler equations (|1.6p as 



povl + Jt{plJ-^),k = Q inf7x(0,T], (1.8a) 

(r;,i;) = (e,uo) in f7 x {i = 0} , (1.8b) 
p'f'^^Q onT , (1.8c) 



with Pq {x) > C dist(a;, F) for x ^ near T. 

1.4. Setting 7 = 2. It should be clear from equations (|1.8p that by introducing new vari- 
ables for both Pq~^ and J^~^, such as the enthalpy for example, we can always return to 
the case that 7 = 2. Henceforth, we seek solutions rj to the following system: 

povl + a'l{plJ-^),k =0 in r! X (0, T] , (1.9a) 

{ri,v) ^ {e,uo) onf7x{i = 0}, (1.9b) 
po==0 onT, (1.9c) 

with po{x) > C dist(a;, F) for x £ Q near F. 
The equation (|1.9ap is equivalent to 

vl + 2A'lipoJ-'),k^0, (1.10) 

and (ll.lOp can be written as 



vl + poaW^:k +Po,3 alJ-'^ = . (1.11) 

Because of the degeneracy caused by po = on F, all three equivalent forms of the com- 
pressible Euler equations are crucially used in our analysis. The equation (|1.9ap is used for 
energy estimates, while (ll.lOp is used for estimates of the vorticity, and (|l.lip is used for 
additional elliptic- type estimates used to recover the bounds for normal derivatives. 

1.5. The reference domain H.. To avoid the use of local coordinate charts necessary for 
arbitrary geometries, for simplicity, we will assume that the initial domain C M"^ at time 
t = is given by 

n^{{xi,x2,x3)eR'' I {xi,x2)eT^, xse (0,1)}, 
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where denotes the 2-torus, which can be thought of as the unit square with periodic 
boundary conditions. This permits the use of one global Cartesian coordinate system. At 
t — Q, the reference vacuum boundary is the top boundary 

r = {x3 = 1} , 

while the bottom boundary {2:3 = 0} is fixed with boundary condition 

m3 ^ on {xa = 0} X [0, T] . 
The moving vacuum boundary is then given by 

T{t)^'n{t){T)^r,{xuX2,l,t). 

1.6. The higher-order energy function. The physical energy /^[^PobP + p^J^'^^dx is 
a conserved quantity, but is far too weak for the purposes of constructing solutions; instead, 
we consider the higher-order energy function 

4 4 
Eit) = J2 Wd^'vml-a + E [\\Pod*-''d?^Drj{t)\\l + W^B^-'^d^^tWo 

a=0 a=0 

+ E llPoa^J-^WIlt^ + 11 curl, vml + WpoB' curl, vml , (1-12) 

where d = (^gf^, gf^) ■ Section [2] explains the notation. 

While this function is not conserved, it is possible to show that supjgjQ E{t) remains 
bounded for sufficiently smooth solutions of (|1.9p , whenever T > is taken sufficiently small; 
the bound depends only on E(0). 

1.7. Main Result. 

Theorem 1.1. Suppose that r]{t) is a smooth solution of m.9\) on a time interval [0,T]. 
Then for < T < Tq taken sufficiently small, the energy function E{t) constructed from the 
solution r](t) satisfies the a priori estimate 

sup E(t) < Mo , 

t6[0,T] 

where Mq and Tq is a function of E(0). 

Of course, our theorem also covers the case that $7 C M'' for d = 1 or 2, and by using a 
collection of coordinate charts, we can allow arbitrary initial domains, as long as the initial 
boundary is of Sobolev class H^-^. We announced Theorem 1 1.1 1 in [4]. 

1.8. History of prior results for the compressible Euler equations with vacuum 
boundary. We are aware of only a handful of previous theorems pertaining to the existence 
of solutions to the compressible and inviscid Euler equations with moving vacuum boundary. 
Makino [I6| considered compactly supported initial data, and treated the compressible Euler 
equations for a gas as being set on M'^ x (0, T]. With his methodology, it is not possible to 
track the location of the vacuum boundary (nor is it necessary); nevertheless, an existence 
theory was developed in this context, by a variable change that permitted the standard 
theory of symmetric hyperbolic systems to be employed. Unfortunately, the constraints on 
the data are too severe to allow for the evolution of the physical vacuum boundary. 

In [llj , Lindblad proved existence and uniqueness for the 3D compressible Euler equations 
modeling a liquid rather than a gas. For a compressible liquid, the density p > is assumed 
to be a positive constant on the moving vacuum boundary T{t) and is thus uniformly 
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bounded below by a positive constant. As such, the compressible liquid provides a uniformly 
hyperbolic, but characteristic, system. Lindblad used Lagrangian variables combined with 
Nash-Moser iteration to construct solutions. More recently, Trakhinin provided an 
alternative proof for the existence of a compressible liquid, employing a solution strategy 
based on symmetric hyperbolic systems combined with Nash-Moser iteration. 

The only existence theory for the physical vacuum singularity that we are aware of can 
be found in the recent paper by Jang and Masmoudi ^6! for the ID compressible gas; we 
refer the interested reader to the introduction in that paper for a nice history of the analysis 
of the ID compressible Euler equations with damping. 

1.9. Generalization of the isentropic gas assumption. The general form of the com- 
pressible Euler equations in three space dimensions are the 5x5 system of conservation 
laws 

p[ut + u ■ Du] + Dp{p) = , (1.13a) 
pt + div{pu) ^ , (1.13b) 
(p(£)t +div(pu(£+pii) = 0, (1.13c) 

where (|1.13ap . (jl.l3b[) and (|1.13cp represent the respective conservation of momentum, mass, 
and total energy. Here, the quantity £ is the sum of contributions from the kinetic energy 
and the internal energy e, i.e.,€ = + e. For a single phase of compressible liquid 
or gas, e becomes a well-defined function of p and p through the theory of thermodynamics, 
e = e{p,p). Other interesting and useful physical quantities, the temperature T{p,p) and 
the entropy S{p,p) are defined through the following consequence of the second law of 
thermodynamics 

T dS — de — dp . 

P 

For ideal gases, the quanities e,T, S have the explicit formulae: 

p T 



e(p,p) 
T{p,p) 



p(7-l) 7-1 
P 
P 

p — p^ , 7 > 1, constant . 



In regions of smoothness, one often uses velocity and a convenient choice of two additional 
variables among the five quantities S, T,p, p, e as independent variables. For the Lagrangian 
formulation, the entropy S plays an important role, as it satisfies the transport equation 

St + {u- D)S = 0, 

and as such, S o r] = So, where So{x) — S{x,0) is the initial entropy function. Thus, by 
replacing / with e^°^ PqJ~'* , our analysis for the isentropic case naturally generalizes to the 
5x5 system of conservation laws. 

2. Notation and Weighted Spaces 

2.1. Differentiation and norms in the open set fl. The reference domain H. is defined 
in Section 11.51 Throughout the paper the symbol D will be used to denote the three- 
dimensional gradient vector 

D - ( — — — 

\dxi ' dx2 ' dx3 
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For integers A: > and a smooth, open domain Q of K"^, we define the Sobolev space 
H''{Q) (il'=(0;]R3)) to be the completion of C°°{Q) {C°°{Q;M.^)) in the norm 



|a|<fc 



1/2 

gai ga2 Qa:, ^^^^ ^ ' 
dXa^ 0X(i2 ^^a3 



for a multi- index a e Z^, with the standard convention that |a| = Oi + 02 + 03. For real 
numbers s > 0, the Sobolev spaces H^{Q) and the norms || • \\s are defined by interpolation. 
We will write H^{Q) instead of H^{Q;R^) for vector-valued functions. In the case that 
s > 3, the above definition also holds for domains O of class H^. 

2.2. Tangent and normal vectors to T. The outward-pointing unit normal vector to F 
is given by 

TV =(0,0,1). 

Similarly, the unit tangent vectors on F are given by 

ri = (1,0,0) and r2 = (0,1,0). 

2.3. Einstein's summation convention. Repeated Latin indices i,j, k,, etc., are summed 
from 1 to 3, and repeated greek indices a, (3, 7, etc., are summed from 1 to 2. For example, 

P „„j pi raf^r'i „ v^2 ^2 dF^ Taff dG* 

^ ■— Z^i=l,3 dxidxi' '« >/3 •— Z^i=l Z^a=l Z^/3=l ax„ dxp- 

2.4. Sobolev spaces on F. For functions u G if'^(F), fc > 0, we set 

1/2 

dx 



\\a\<k-^^ 



u(x) 



dX(xi dXa2 



for a multi-index a e Z^. For real s > 0, the Hilbert space H^{T) and the boundary norm 
I • Is is defined by interpolation. The negative-order Sobolev spaces H~^{r) are defined via 
duality: for real s > 0, 

2.5. Notation for derivatives and norms. Throughout the paper, we will use the fol- 
lowing notation: 

^ n . , (Odd 

D = three-dimensional gradient vector = — — , — — , — — 

\ OXl OX2 OXz 

B = two-dimensional gradient vector or horizontal derivative = I -- — , — — ) , 



9a; 1 8x2 



II • lis = H^{d) interior norm, 
I • Is = H^{r) boundary norm. 

The fcth partial derivative of F will be denoted by F,/. = 
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2.6. The embedding of a weighted Sobolev space. Using d to denote the distance 
function to the boundary T, and letting p = 1 or 2, the weighted Sobolev space iJ^p(ri), 
with norm given by J^d{x)P{\F{x)\'^ + \DF{x)\'^) dx for any F G H^p{Q), satisfies the 
following embedding: 

so that there is a constant C > depending only on Q and p such that 

\\F\\l_^^, < C jj{xY{\F{x)\^ + \DF{x)\') dx. (2.1) 

See, for example, Section 8.8 in Kufner 

3. The Lagrangian vorticity 
We make use of the permutation symbol 

{1, even permutation of {1, 2, 3}, 
— 1, odd permutation of {1, 2, 3}, 
0, otherwise, 

and the basic identity regarding the ith component of the curl of a vector field u: 

(curlu)i = £ijkU^ ■ 

The chain rule shows that 

(curlu(?7))i = curL^u e.ijkA^^v'' ,s , 

the right-hand side defining the Lagrangian curl operator curl^. Taking the Lagrangian curl 
of p.lOp yields the Lagrangian vorticity equation 

ekjiAjvl,s = , or curl^ = . (3.1) 
4. Properties of the determinant J, cofactor matrix a, unit normal n, and a 

POLYNOMIAL-TYPE INEQUALITY 

4.1. Differentiating the Jacobian determinant. The following identities will be useful 
to us: 

dJ = a'Ld— — (horizontal differentiation ) , (4-1) 

ox" 

dtJ — a^— — (time differentiation using v — ijt) . (4.2) 
ox^ 

4.2. Differentiating the cofactor matrix. Using (|4.ip and (|4.2p and the fact that a — 
J A, we find that 

da!l = d-^J^^[ala^ — a^a^] (horizontal differentiation) , (4.3) 
dto-i = -g—^J^^[a'^a^ ~ afa'^] (time differentiation using w = jyt) . (4.4) 

4.3. The Piola identity. It is a fact that the columns of every cofactor matrix are divergence- 
free and satisfy 

a-,fc = 0. (4.5) 

The identity (|4.5p will play a vital role in our energy estimates. (Note that we use the 
notation cofactor for what is commonly termed the adjugate matrix, or the transpose of the 
cofactor.) 
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4.4. Geometric identities. The vectors r],a for a = 1,2 span the tangent plane to the 
surface F in M'^, and 

V,l V,2 , '7,1X77,2 

Ti :— , T2 := -. r , and n :— 



\v,i I ' ' \V,2 1 ' ' ha XV,2 I 

are the unit tangent and normal vectors, respectively, to F. 

Let gaj3 = ri,a -rj^p denote the induced metric on the surface F; then det <? — \r],i xri,2 ^ 
so that 

Vgn := 77,1 X77,2 , 



where we will use the notation y/g to mean y/detg. 
By definition of the cofactor matrix. 



'?^lf7^2-?7^l^7^2 

V^,lV'^,2-V^,2V'^,l 



, and |a^|. (4.6) 



It follows that 



n = a^/V9- (4.7) 



We will often make use of the following differentiation formulas for the unit normal and 
tangent vectors: 

fit = -g'^'^{v,i3 -n)!],^ , 

where g^^ denote the inverse of the metric g^p. Note that the right-hand sides these identities 
are tangent vectors to the embedded surface. 

4.5. A polynomial-type inequality. For a constant Mq > 0, suppose that f{t) > 0, 
1 1-^ f{t) is continuous, and 

fit)<Mo + CtPif{t)), (4.8) 

where P denotes a polynomial function, and C is a generic constant. Then for t taken 
sufficiently small, we have the bound 

fit) < 2Mo . 

This type of inequality, which we introduced in can be viewed as a generalization of 
standard nonlinear Gronwall inequalities. 

With E{t) defined by (|1.12p . we will show that supj^jQ i?(i) satisfies the inequality 



5. Trace estimates and the Hodge decomposition elliptic estimates 



The normal trace theorem which states that the existence of the normal trace of a velocity 
field w G L^(il) relies on the regularity of divw (see, for example, [E]). If divw G H^ifl)', 
then w ■ N, the normal trace, exists in H^'-^-^{T) so that 



||w • iV||^-o.5(r) < C ||w||i2(o) + ||divw||^i(j^), 



(5.1) 



for some constant C independent of w. In addition to the normal trace theorem, we have 
the following. 



FREE-BOUNDARY 3D COMPRESSIBLE EULER EQUATIONS IN VACUUM 



9 



Lemma 5.1. Let w G L^(r2) so that curlw G H^{fiy , and let Ti, T2 denote the unit tangent 
vectors on T, so that any vector field u on T can be uniquely written as u^Ta- Then 



hllL(o) + l|curlz«||^i(n)'J > a =1,2 (5.2) 
for some constant C independent of w. 

See [1] for the proof. Combining (|5.ip and (|5.2p . 



\w\\h-«'-(t) 



< C 



\w\\L2{n) + ||divw||Hi(o)' + l|curlu;||Hi(o)' (5.3) 



for some constant C independent of w. 

The construction of our higher-order energy function is based on the following Hodge- type 
elliptic estimate: 

Proposition 5.2. For an H'' domainft, r > 3, ifF G L'^{n;R^) withcmlF G H^-^iniM.^), 
divF G 7J'*~^(ri), and F ■ N\r G if'*"^(r) forl<s<r, then there exists a constant C > 
depending only on fl such that 

\\F\\s < C (||F||o + II curl^^||,_i + II divF||,_i + \dF ■ TV^^s) , 

Ili^lU < C (\\F\\o + II curlF||,_i + II divF||,_i + \dF ■ T^l^.a) , 

where N denotes the outward unit-normal to T, and Ta are tangent vectors for a = 1,2. 

The first estimate is well-known and follows from the identity — AF = curlcurlF— Z3divi^; 
a convenient reference is Taylor [17| . The second estimate follows from the first using the 
same geometric identities on the boundary. 

6. The a priori estimates 

Since the degeneracy of the initial density is only in the normal (or vertical) direction to 
the vacuum boundary, and hence there is a constant C > such that |9po(2;)| < Cpo{x), we 
may assume without loss of generality that pQ = po{x3) and po,3 (2^3) = 1 for very small. 
In fact, it is convenient to suppose that 

^0(2^3) = 1 - a;3 , 

although any sufficiently smooth function which vanishes on T and is bounded from below 
by a constant multiple of the distance function near F would suffice. 

6.1. Curl Estimates. Following Lemma 10.1 in [3], we obtain the following estimates. 

Proposition 6.1. For all t G (0,T), 
3 4 
Y,\\cni:\d^''vml-a+T.\\Po9'-'cmld^'r^{t)\\l < Mo + CTP{ sup E{t)) . (6.1) 

Proof. From (|3.ip . (curl^w)j' = EkjiAtjV'^^s B{A,Dv), where B is quadratic in its argu- 
ments; hence. 



curl,, v{t) = curl Wo + / B{A{t'), Dv{t'))dt' , (6.2) 
Ja 

and computing the gradient of this relation yields 

curl^ Dv{t) = curlDuo - e.j.D^^u*,^ + [ DB{A(t'),Dv{t'))dt' . 

Jo 
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Applying the fundamental theorem of calculus once again, shows that 

pt pt i*t 

curl„DT]{t) =tcuT\Duo+e.ji [At]Dr]\s -DA'^v'^sW + / DB{A{t"), Dv{t"))dt"dt' , 

Jo Jo Jo 

and finally that 

curlL»77(t) = i curl Duo - J At]{t')dt' Drf 

+ e.j,f[At'^Dif,s~DA'^v\s]dt'+f f DB{A{t"),Dv{t"))dt" dt' . (6.3) 
Jo Jo Jo 

To obtain an estimate for || curl?7(i)|||, we let D'^ act on ([131). With dtA'^ = -Afv^^A^. 
and DA^j = —AfDrf^p A^, we see that the first three terms on the right-hand side of (|6.3p are 
bounded by Mq + CTP(sup(gjo t] where we remind the reader that A/q = P{E{Q)) 

is a polynomial function of the E at time t — 0. Since 

DB{A, Dv) = -£uJ^[Dv\, Atv\^ A^ + v\, AfDv^^ A^^ + v\, v',p DiAfAP)], 
the highest-order term arising from the action of on DB{A, Dv) is written as 

'£kJ^ f f[D^v\sAtv\pAP + v\,AtD^v\pAP]dt"dt' . 
Jo Jo 

Both summands in the integrand scale like D^vDvAA. The precise structure of this sum- 
mand is not very important; rather, the derivative count is the focus. Integrating by parts 
in time. 



PL PZ PZ PZ PZ 

/ / D^vDvAAdt"dt' = ~ / D^T]{DvAA)tdt"dt'+ D^?]DvAAdt' 
Jo Jo Jo Jo Jo 

lows that 

pt pt' 

/ / D'^B{A{t"),Dv{t"))dt"dt'\\l<CTP{ sup E{t)), 

Jo Jo t6[0,Tl 



from which it follows that 
pt pt' 



and hence 

sup ||curh/(t)||^ < Afo + CTP( sup E{t)) . 
te[a,T] te[o,T] 

Next, we show that 

II curlut(<)||^ < Mo + CTP{ sup E{t)) . (6.4) 

te[o,T] 

From (jXTI) . 



cnv\vt=e.j., / At'^{t')dt' V, 
Jo 

Since H'^ (fi) is a multiplicative algebra, we can directly estimate the ^ (ri)-norm of curl Vt 
to prove that (|6.4p holds. The estimates for curltittt(i) in H^[n) and curl9f'i;(t) in L^(r2) 
follow the same argument. 

The weighted estimates follow from similar reasoning. We first show that 

||po^^curl7?(<)||2 < Mo + CTP( sup E{t)) . (6.5) 

te[o,T] 
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mrlvit) = £jkiv\s ( At'{t')dt' + cnrluo+ [ B{A{t'), Dv{t'))dt' , 
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and integrate in time to find that 



curl ry(i) = t curl Mo + / ejk^v\sf Atj{t")dt"dt' + [ [ B{A{t"), Dv{t"))dt"dt' 
Jo Jo Jo Jo 



It follows that 
Pod^ cur\ri{t) — tpoB'^ cuAuq 

rt rt' 

ekJ^At'po^^v\,dt"dt' + 



JO 

t 



t rt' 



SkJ^Po^ Aty,sdt"dt' 



+ I SjkiPod'^v' 



AtUt")dt"dt' + / eJk^v\s / pod^AtUt")dt"dt' + ^2, (6.6) 



where 9^2 denotes terms which are lower-order in the derivative count; in particular the 
terms with the highest derivative count in scale like pd^Dv or pd'^rj, and hence satisfy 
the inequality |19^2(i)||o ^ -^o + C7'^'(suptg[o,T] ^(0)- We focus on the first integral on the 
right-hand side of (|6.6p : integrating by parts in time, we find that 



t rt' 



£kjiAtjPod'^v\s dt"dt' 



10 Jo 
and hence 



t rt' 



Jo 



£kjiAu]pod^rf,sdt"dt'^ j EkjiAt^pod^if ,s dt' 



t rt' 



Jo 



£k3^At'pod^v\sdt"dt' 



<Mo + CTP{ sup E{t)). 

te[o.T] 



The other time integrals in (|6.6[) can be estimated in the same fashion, which proves that 
(|6.5p holds. The weighted estimates for the curl of Vt, Vut and d^v are obtained similarly. □ 

6.2. Energy estimates. We assume that we have smooth solutions 77 on a time interval 
[0,T], and that for all such solutions, the time T > is taken sufficiently small so that for 

t e [0, T] 

Mt)\\i,<2\h+\' + l, (6.7) 
\\dMml-a/2<ndMml-a/2 + '^ for a = 0,1,..., 6. 

The right-hand sides appearing in these inequalities shall be denoted by a generic constant C 
in the estimates appearing below. Once we establish our a priori bounds, can indeed verify 
that our solution adhere to the assumptions (|6.7p by means of the fundamental theorem of 
calculus. 



6.2.1. The structure of the estimates. Due to the degeneracy of the initial density function 
Poi one time derivative scales like one-half of a space derivative. The energy estimates for 
the time and tangential derivatives are obtained by first studying the ^^-differentiated Euler 
equations, then the 9'^9t^-differentiated Euler equations, and so on, until we reach the d^df- 
differentiated Euler equations. The estimates for the normal derivatives are then found 
using elliptic- type estimates. The Sobolev embedding theorem requires that we use H'^{Vl) 
as the minimal regularity of 77 (i). 
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6.2.2. The d'^ -problem. 

Proposition 6.2. For S > and letting the constant Mq depend on l/S, 

sup ( / po{x)\d'^v{x,t)\^dx+ I pl{x)\d^Dri{x,t)\^dx 
te[o,T] \JQ. Jn 

<Mq + S sup E{t) + CTP{ sup E{t)). (6.8) 
tG[0,T] te[o,T] 

Proof. Letting d"^ act on pqvI + a^(pgj^^),fe = 0, and taking the L^(f7)-inner product with 
d'^v'', we obtain 



I P^\dWdx+ I d^a'l{plJ-^),kd\'dx+ I a'y{pld\j-^),kd%'dx 

3 



Integrating the first term from to i G (0, T] produces the first term on the left-hand side 
of (Oil . 

We define the foUowing three integrals 



Ii = / d''at{ptJ-'),kd\'dx 
Jn 

l2= f a'y{pld\j-''),ud\'dx 
Jn 

3 



1=1 



The last integral introduces our notation TZ for the remainder^ which throughout the paper 
will consist of integrals of lower-order terms which can, via elementary inequalities together 
with our assumptions (|6.7p . easily be shown to satisfy the following estimate: 

/ n{t)dt<Mo + S sup E[t) + CTP{ sup E{t)). (6.9) 
-'0 te[o,T] te[o,T] 



The sum of \Ii{t) +l2{t)]dt together with the estimates for curl 77 given by Proposition 
16.11 will provide the remaining energy contribution j^^p^{x,t)\d'^Dri\^dx plus error terms 
which have the same bound as TZ. 

Analysis of jj" TZdt. We integrate by parts with respect to x^ and then with respect to 
the time derivative dt, and use (|4.5p to obtain that 



7^ = -Vc^ /" ( d^-'a'l pIB^J-^ d\\k dxdt 
1=1 "'0 "'^ 

3 f r 

= Vq / / po{d^-'a'^d'j-^)^pod^7j\kdxdt-y2ci / pod^-'a'^d'j-^podW,kdx 
1=1 -^0 Jn Jn 



Notice that when ^ = 3, the integrand in the spacetime integral on the right-hand side 
scales like £ [dDrj pod^dtJ~^ + dDv pod^J~^] pod^Drj where £ denotes an L°°(ri) function. 
Since \\podfJ~^{t)\\-^ is contained in the energy function E{t) and since dDr]{t) G L°°{Q), 
the first summand is estimated using an L°°-L'^-L'^ Holder's inequality, while for the second 
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summand, we use that Upo^/^^lOlll is contained in E{t) together with an L'^-L'^-L^ Holder's 
inequahty. 

When / = 1, the integrand in the spacetime integral on the right-hand side scales like 
£ [dDri poB'^ati + dDv poB^a^] p^d'^'q^.k. Since \\pad'^ Dvt{t)\\Q is contained in the energy 
function E{t) and since dD-q e L°°(f7), the first summand is estimated using an L°°-L^-L^ 
Holder's inequality. We write the second summand as 

BDvpod^af podS\.p+dDvpod\f poB^v^s ■ 

We estimate 

/ BDvpqB^4 PoB*f]\0 dxdt = - ( ( [BDvpoB^a^,p pqB^t]' + BDv,i3 poB^a^ pQB^f]']dxdt 
Jn Jo Jn 

<C [ {\\BDv{t)\\Ls^n)\\poB*a{t)\\o WpoB^vmL^n) 
Jo 

+ \\B^Dv{t)\\Ls^n)\\poBSimLein)\\B^a\\o)dt 

<C [ {\\BDv{t)\\Ho..^n)\\poB^a{t)\\o WpoB^vimi 
Jo 

+ \\B^Dv{t)\\Ho.^n)\\poB^vmi \\B^a\\o)dt 

<C f {\\vmH^-Hn)\\poB^Dr^{t)\\l + Mt)\\H--Hn)\\poB^D7^moMt)h 
Jo 

+ \\v{t)\\H^.^n)\\vml)dt, (6.10) 

where we have used Holder's inequality, followed by the Sobolev embeddings 

H°-^{n) ^ L^{n) and H\n) ^ L^ifl) . 

We also rely on the interpolation estimate 
2 
ll"llL2(o,T;ff3-5(n)) ^ C'(||u(t)||3||r7||4) ^ + C\\vt\\ L^(^o^T:min))\\v\\ L^iO,T;m{(l)) 

<Mo + 6 sup Ut)\\l + CT sup (Mt)\\l + \\vtml). (6.11) 
te[o,T] te[o,T]^ ^ 

where the last inequality follows from Young's and Jensen's inequalities. Using this together 
with the Cauchy-Schwarz inequality, (|6.10p is bounded by CrP(suptg[o t] ^(0)- Next, 
since (14. 6p shows that each component of is quadratic in Brj, we see that the same 
analysis shows the spacetime integral of BDv poB^a^ p^B'^if has the same bound, and so 
we have estimated the case ^ = 1. 

For the case that 1 = 2, the integrand in the spacetime integral on the right-hand side 
of the expression for TZ scales like i B^ Drj B^ Dv poB'^Drj, so that an — — Holder's 
inequality, followed by the same analysis as for the case I = 1 provides the same bound as 
for the case ^ = 1. 

To deal with the space integral on the right-hand side of the expression for TZ, the integral 
at time i = is equal to zero since r]{x, 0) = x, whereas the integral evaluated at t = T is 
written, using the fundamental theorem of calculus, as 

3 37^ 

-Vq / poB^~'a'yB'j~^poB^Tj\kdx =-Vq /po / iB^~'4B'j-%poB^Tj\k{T)dx 
i^-^ Jn *='^ Jn Jo 

which can be estimated in the identical fashion as the corresponding spacetime integral. As 
such, we have shown that TZ has the claimed bound 
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Analysis of the integral Ii. Because po = on F = {xa — 1}, we integrate by parts to 
find that 

Ii = - I plJ-^d'^a'yd'^v\kdx+ [ pp-^ d'^afd^v' dxidx2 
Jn J{x3=o} 

= - I plJ-^B^alBWkdx, 
Jn 

since on the fixed boundary {a;3 = 0}, rj^ = 0:3 so that according to (|4.6p . the components 
af = and = on {2:3 = 0}, and = on {xs = 0}, so that B^a^d'^'v^ = on{x3 = 0}. 

To estimate Ii, we use the formula (|4.3p for horizontally differentiating the cofactor 
matrix: 

Ii = / po^J-^BSf,s Ka^-<af] 5V,fc dx + n, 
Jn 

where the remainder TZ satisfies ()6.9p . We decompose the highest-order term in Ti as the 
sum of the following two integrals: 

Iia= [ (aV,saD(5V,fea,^)dx, 
Jn 

In = - f (9V,. <)(^V,fea,'=)dx. 

Jn 

Since u = 774, Ii^ is an exact derivative modulo an antisymmetric commutation with respect 
to the free indices i and r; namely, 

d^7^\s ald\\k a': = dS\s KBWk 4 + {d*v\s at - B^ K)B\\k 4 . (6.12) 
Using the notation 

[DrjFjl = alF\s for any vector field F , 

BS\sa:.B\\kal = —WW" - iaV,. ~dS\k {<a1)t , (6.13) 

so the first term on the right-hand side of (|6.12p produces an exact derivative in time. 
For the second term on the right-hand side of (|6.12p . note the identity 

at - BSf,s at)B\\k = -Jh^JkBS^r A] e^mnB\'',s . (6.14) 

We have used the permutation symbol e to encode the anti-symmetry in this relation, and 
the basic fact that the trace of the product of symmetric and antisymmetric matrices is 
equal to zero. 

Recalling our notation [curh,F]* = £^7=^'=,^ AJ, (l6l4)l can be written as 

at - B^rf,s at)B\\k = - curl^ 9^/ • curl^ B\ , (6.15) 
which can also be written as an exact derivative in time: 

curl, BS ■ curl, B^v = curl, B^l' - iA-^A'^)t + B^^r BS\s {A]Al)t . 

(6.16) 

The terms in (j6.13p and (j6.16p which are not the exact time derivatives are quadratic in 
Pod'^Drj with coefficients in L°°{[0,T] x fi); denoting the integral over fl of such terms by 

^^a = ~ [ Po'J-'\D^B'vfdx-~ I p^J-^\cuT\BS\^dx + Q^^Q.j,^+Tl, 
where Jq \Qpod'^Dr,\ dt < C T P{snpte[o,T] E{t)), and TZ satisfies (jej]) . 



FREE-BOUNDARY 3D COMPRESSIBLE EULER EQUATIONS IN VACUUM 



15 



With the notation div^ F — A^F'^jj, the differentiation formula (14. 1|) shows that In, can 
be written as 



1- d f 



It follows that 



■ f po'(J"^|-D9\|2- J-i|curl„a\p-J-i|div^9\p)d.T + 7e, 
Jn 



2di 



where we have used the fundamental theorem of calculus for the second equality on the term 

Analysis of the integral T2. Integration by parts once again yields 

2:2 = - / pldU-^4d%\kdx. 
Jn 

Since d'^J^^ ~ —2J^^d*J plus lower-order terms, which have at most three horizontal 
derivatives acting on J. For such lower-order terms, we integrate by parts with respect to 
dt, and estimate the resulting integrals in the same manner as we estimated the remainder 
term TZ, and obtain the same bound. 
Thus, 



l2 = 2 [ plJ'^ald'^r]',r 4d*v\kdx + n 
Jn 



/n 

= H Pp'^'aldS^r a^dS\kdx- [ pl{J-^ala\)t B^r^' ,r ,k dx + TZ 
Jn Jn 

Given our identities for differentiating a and J, the Sobolev embedding theorem together 
with our assumptions (j6.7p and the Cauchy-Schwarz inequality show that 

/ / pl[J'^ala\)tB^T]\rBS\kdxdt<CT sup E[t)- 
Jo Jn te[o,T] 

consequently, we can write 

rt 

o2 7-3„ro4^s „fea4„ 



pp-^ald^T^',^ a^d^r]\udx^ Ma+ [l2{t') + n{t')]dt' . (6.17) 

! Jo 

On the other hand, 

plJ-^a:BS',r a'ldW,kdx 

= f pIJ~'^{B^ diY 7] + BS',r f atldt') {B'^diYri + B'^T^\k I at^dt')dx 
Jn Jo Jo 

= f plJ-'^\B^divr]\^dx + 2 f pp^^B* div rj 8^' ,r f at^dt' dx 
Jn Jn Jo 

+ I pIJ-^ BS',r I atldt' BWk I at\dt' dx (6.18) 
Jo Jo Jo 

Yet another application of the Sobolev embedding theorem together with our assumptions 
(|6.7p and the Cauchy-Schwarz inequality shows that the second and third integrals on the 
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right-hand side are bounded by Mq + CT suptg[Q so that combining (I6.17P and 

(Hm), we find that 



pp-^ld* diYr]\'^dx ^ Mo + / [J2it') + n{t')]dt' 



(6.19) 



Summing inequalities. Summing ()6.19|) with Ti yields 



plJ-^\d^Dri\^dx+ / plJ-^ld'^divTjl^dx- / plJ-^\d'^ cmlT]\^dx 



<Mo + 5 sup E{t) + CTP{ sup E{t)) . 
te[o,T] te[o,T] 



sup - 

te[o,T] ^ 



Adding to this, the inequahty (16. and possible readjusting our constants, we obtain the 
desired result, and complete the proof of the proposition. □ 

Since 77 • = 77" for a = 1, 2, we have the following 

Corollary 6.3. For a = 1, 2, 

sup |?7"|^ 5 < Mo + CTF( sup E{t)) . 
tG[0,T] te[o,T] 

Proof. The weighted embedding estimate (|2.ip shows that 



Now 



\\BS\\l<C pl{\dS\^ + \d'Dr^\')dx. 



te[o,T] Jn 



sup / pQ|9^?7pdx = sup 



tG[0,T] Jo 







2 


/ pI 


[ d^vdt' 


dx<T^ sup 


in 


Jo 


*G[0,T] 



It follows from Proposition 16.21 that 



sup \\d^v\\l< Mo + CTP{ sup E{t)). 



te[o,T] 



te[o,T] 



According to our curl estimates (|6.ip . supj^ro x] |1 curl?7||| < Mo + C T P{supf.^iQ j,, E{t)), 



from which it follows that 



sup ||a''curl?7||^i(f^), < Afo + CTP( sup E{t)) , 



tG[0,T] 



te[o,T] 



since 9 is a horizontal derivative, and integration by parts with respect to d does not produce 
any boundary contributions. From the tangential trace inequality (|5.2p . we find that 



sup |(9Vl-i/2 < *^o + CTP( sup E{t)), 
te[o,T] *e[o,T] 

from which the assertion of the corollary follows. 
6.3. The 9f -problem. 

Proposition 6.4. For S > and letting the constant Mq depend on l/S, 
po\d^v{x,t)\'^dx+ / pl{x,t)\djDv{x,t)\'^dx 



□ 



sup 

tG[0,T] \Jn 



<Ma + 5 sup E{t) +CTP{ sup E{t)) . (6.20) 
te[o,T] te[o,T] 
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Proof. Letting df act on povl + a'[{plJ^^),k — 0, and taking the L^(fi)-inner product with 
dfv"^, we obtain 



2di 



Po\dfv\'dx+ / dX{pU-'),kdydx+ / aUpidfJ^'Ud^dx 
n Jn Jn 



j^ci f dr'aUpld\j-^).kdx. 
1=1 



Integrating the first term from to t G (0, T] produces the first term on the left-hand side 
of (E^Ol)- 

We define the following three integrals 



Xi= / d^,at{ptJ-'),udtv^dx 
Jn 

l2= f a'y{pld^J~^),kd!v'dx 
Jn 

n = Y^ci[ df-'a1 {pld\.r^),u d!v^ dx . 

, -, Jn 



1=1 

The sum of J^[Ji{t) +J2it)]dt together with the curl estimates given by Proposition 16.11 
will provide the remaining energy contribution Po{x, t)\dj Dv\'^ dx plus error terms which 
have the same bound as TZ. 

Analysis of J^^ TZdt. We integrate by parts with respect to Xk and then with respect to 
the time derivative dt to obtain that 



7 X 

n = -Y,ci f f d^-^a'i; pldlJ-^ dfv\k dxdt 
, Jo Jn 



1=1 

7 „T 



J2ci f f po {dt'a'ldlJ-'),PodJv\k dxdt - ^ Q / 
1=1 ■^^ -^^ 1=1 



Pod!-'a'id'j-^podJv\kdx 



T 



Notice that when / — 7, the integrand in the spacetime integral on the right-hand side 
scales like £ [Dvt podfDv + Dv podjDv] podjDv where i denotes an L°°{n) function. Since 
\\podjDv{t)\\g is contained in the energy function E{t), Dvt{t) is bounded in L°°(r2), and 
since we can write podfDv{t) = podfDv{0) + podj Dv{t')dt' , the first and second sum- 
mands are both estimated using an L°°-L^-L^ Holder's inequality. 

The case Z = 6 is estimated exactly the same way as the case / = 3 in the proof of 
Proposition 16.21 For the case / = 5, the integrand in the spacetime integral scales like 
^[DvttPQ^^ J~'^ + DvutPoDvtttt]PodJ Dv. Both summands can be estimated using an L^-L^- 
Holder's inequality. The case I = 4 is treated as the case / = 5. The case Z = 3 is also 
treated in the same way as I = 5. The case I = 2 is estimated exactly the same way as the 
case Z = 1 in the proof of Proposition 16. 21 The case Z = 1 is treated in the same way as the 
case I = 7. 

To deal with the space integral on the right-hand side of the expression for TZ, the integral 
at time t = is bounded by Mq, whereas the integral evaluated at f = T is written, using 
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the fundamental theorem of calculus, as 

7 



/ podl~'a1dU-^podlv\kdx =y^ci I p„dt'a^{0)dlJ-'{0)podJv\kiT)dx 
1=1 •'^ 1=1 

j^-^ Jo Jo 

The first integral on the right-hand side is estimated using Young's inequality, and is 
bounded by Mq + 5 sup^gjo t] -^(0' while the second integral can be estimated in the iden- 
tical fashion as the corresponding spacetime integral. As such, we have shown that TZ has 
the claimed bound (16.91). 



Analysis of the integral 2i. As to the term Xi, using the identity (|4.4|) . the same 
computation as for the ^^-differentiated problem shows that 

pMv\,AI) {dlv\^A':) = ]^^^\p,D^dlv{t)\^ l||p„curl,9j«(0P 

+ \p^dlv\rdlv\, {A^^At^mjt - 5161] , 



and 



-plidy,, AD {dtv\uA':) = -iA|podiv,5M2 + ]^pldy,,dlv\u (KA% 



and hence 



= / po^ [j~'^\D,p'lv\^ - J-^\cm\^d'lv\^ - J-^\diY^d'lv\^)dx + 'R. 
2 dt Jfi 

Analysis of the integral T2. Integration by parts once again yields 

^2 = - [ pld^J-^a'idfv\kdx. 
Jn 

Since dfj~'^ = —2J~'^dfj plus lower-order terms, which have at most seven time derivatives 
on J, and can be estimated in the same fashion as the remainder term TZ above. 
We see that 

J2 = 2 /" pp-^aldlv^.r a'^dfv\kdx + n 
Jn 



d_ 

dt 



' pp-^aldjv^r a^djv\kdx- f pliJ~''ala'^)t djv%r djv\k dx + TZ 
n Jn 

4 



Following our analysis of the term X2 in the 9 -problem, we see that 

plJ-^aldJv',r a'ldJv\kdx^Mo+ f [I2 {f )+ n{t')]dt' . (6.21) 
n Jo 
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On the other hand, 

plJ-^{dJdivv + dJv',r I a/Jt') {dj divv + djv\k I atfdt')dx 



f plJ~^\dJdivvfdx + 2 f plJ~^dJ divv djv' ,r f at^dt' dx 
Jn Jn Jo 

pIJ-^ djv',r f atldt' djv\k f at'^dt' dx (6.22) 



n 



Yet another apphcation of the Sobolev embedding theorem together with our assumptions 
()6.7p and the Cauchy-Schwarz inequahty shows that the second and third integrals on the 
right-hand side are bounded by Mq + C supjgjg.T] so that summing (I6.2ip and (I6.22|) 

shows that 

/ plJ-^\dJ divv fdx ^ Mo + [ [l2{t')+n.{t')]dt' . (6.23) 
Jn Jo 

Summing inequalities. Summing (|6.23p with Ti yields 



1 r 
sup t; 

te[o,T] ^ 



[ plJ^^\djDv\'^dx + plJ-^\dJ divv\'^dx - [ pp^^ldj cnTlv\^dx 
Jn Jn 

< Mo + 6 sup E{t) + CTP{ sup E{t)) . 
te[o,T] «e[o,T] 



Adding the curl estimate (|6.ip . readjusting our constants, we obtain the desired result, and 
complete the proof of the proposition. □ 

6.4. The dfd^, df^^i '^t^ problems. Since we have provided detailed proofs of the 
energy estimates for the two end-point cases of all space derivatives, the B"^ problem, and 
all time derivatives, the df problem, we have covered all of the estimation strategies for all 
possible error terms in the three remaining intermediated problems; meanwhile, the energy 
contributions for the three intermediate are found in the identical fashion as for the 9* and 
df problems. As such we have the additional estimate 

Proposition 6.5. For S > and letting the constant Mo depend on l/S, for a — 1,2, 

sup ^[■^'"'^"Wla.s-a + WVf^d^-" df'vml + \\pod'-'^ d^^Drj{t)\\l 

<Mo + S sup E{t) + CTP{ sup E{t)) 
te[o,T] te[o,T] 

6.5. Additional elliptic-type estimates for normal derivatives. Our energy estimates 
provide a priori control of horizontal and time derivatives of r/; it remains to gain a priori 
control of the normal (or vertical) derivatives of rj. This is accomplished via a bootstrapping 
procedure relying on having djv{t) bounded in L^(f2). 

Proposition 6.6. For t G [Q,T], dfv{t) e H^{n), podfj-'^it) e H\n) and 

sup {\\d^v{t)\\l + \\podU-\t)\\l) <Mo + 6 sup E{t) + CTP{ sup E{t)) . 
te[o,T] te[o,T] te[o,T] 
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Proof. We write (|1.9ap as vl + 2A^{poJ ^),k ~ 0, which we rewrite as 

vl + poa- -2af = . (6.24) 

We have used the fact that po,/3 — for /3 — 1, 2, and po,3 = —1- Letting act on equation 
((04)) . we have that 

5 

a=l 

According to Proposition 16.41 and 16.51 

sup (\\dJv{t)\\l + \\podDd^v{t)\\l) < Mo + S sup E{t) + CTP{ sup £;(i)) , 
te[o,T]^ ^ te[o,T] te[o,T] 

and since (|4.6p shows that is quadratic in c??/, we see that for ah t e [0,T], 

\\[poa^^dfj~^,3-2afd^J-^]{t)\\l < Mq + 5 sup E{t) + CTP{ sup £:(t)) . 

te[o,T] te[o,T] 

It foUows that 

||pol«:'|5fJ-^3(<)ll§ + 4|||a3|9,V-2(t)||2-4 / poW^l'dfj-'dU-' ,3 dx 



<Mo + S sup i;(t) + CTP{ sup £;(i)) . 
tG[o,r] te[o,T] 

We assume that our solution is sufficiently smooth so that dfj^^{t) e H'^{n), and in 
particular that J^^)^],3 is well-defined and integrable. As such, we writ^H 

-4 / polar'I'a.V-^Sfj-^^g dx^-2\\\^\dU-Ht)\\l + ^ I Po(|a'n,3 (aej-^)^^^ 
Jo Jo 

+ 4 /" |9f J"2|2^2;ida;2, 
"'{3:3=0} 

so that together with our previous inequality, 

\\p,du~\3 {t)\\i + \\du-\mi 

<Mo + S sup E{t) + CTP{ sup E{t)) + C [ po\d'iJ-^\^dx. 

te[a,T] te[a,T] Jo 

Since poddfJ^^{t) is already estimated by Proposition 16.51 then 

\\pod^j-'mi + \\d!j~'mi 

<Mn + S sup E{t) + CTP{ sup E{t)) + C [ po\d^ J^^f dx . 
te[o,T] te[o,T] Jo 



^Jang & Masmoudi ^ have counterexamples to the obtained inequahty when is not sufficiently 

smooth. It is important that the function has greater regularity than the desired a priori estimate 

indicates, and in particular, as we noted, J~^)^],3 must be well-defined and integrable. 



FREE-BOUNDARY 3D COMPRESSIBLE EULER EQUATIONS IN VACUUM 



21 



We use Young's inequality and the fundamental theorem of calculus (with respect to t) 
for the last integral to find that for 6 > 

C [ podfj-^d^J-^dx<S\\dfj-^{t)\\l + Cs\\pod^^Dvit)\\l+Mo + CTP{ sup E{t)) 
Jq+ ie[0,T] 

<S\\dfj-^{t)\\l + Mo + CTP{ sup Eit)), 

te[a,T] 

where we have used the fact that \\podJ Dv{t)\\Q is contained in the energy function E{t). 
By once again readjusting the constants, we see that on [0,T] 

\\podfj-^it)\\l + \\d^J-^{t)\\l <Mo + S sup E{t)+CTP{ sup E{t)) . (6.25) 

te[o,T] te[o,T] 

With Jt — ajv^jj, we see that 

4 

a=l 

SO that using (|6.25p together with the fundamental theorem of calculus the estimate for the 
last two terms on the right-hand side, we see that 



ald'lv'^j (t) < Mo + S sup E{t) + CTP{ sup E{t)) , 



te[o,T] 



te[o,T] 



from which it follows that 



|divafw(t)|| <A/o + (5 sup E{t) + CTP{ sup E{t)) . 



te[o,T] 



te[o,T] 



According to Proposition 16. li || curl(?fw(i)||o < Mq + CTF(suptg[Q ^]i?(t)) and with the 
bound on d^v" given by Proposition 16.51 Proposition 15.21 provides the estimate 



|afi;(t)||:; < Mq + 6 sup E{t) + CTP{ sup E{t)) . 

tG[0,T] te[o,T] 



□ 



Having a good bound for dfv{t) in H^{il.) we proceed with our bootstrapping. We let 
act on equation (|6.24p . so that 

poaldfj-^, -2afdtJ-' = ~d!v^ - podfia^J-^p ) - {dfal)[-2J-' + po J-^3 ] 

3 



with the right-hand side bounded in H^{n) by Afo + <5suptg[Q rp^E{t) + C T P{s\ip^^^Q yji?(t)). 
Using the argument just given, we conclude that 

sup {\\vtttml + \\podtJ-Ht)\\l) <Mo + S sup E{t) + CTPi sup E{t)) . 
te[o,T] te[o,T] te[o,T] 



Next, we let act on equation (|6.24p . so that 



+ 2{dtal)dt[~2J-^ + poJ-\i], 



22 



D. COUTAND, H. LINDBLAD, AND S. SHKOLLER 



with the right-hand side bounded in H'^{Q) by A/o + <5suptg[Q + CTP(supjg[Q 

We conclude that 

sup {\\vt{t)\\l + \\pod^J-'^{t)\\l) <Mo + d sup E{t)+CTP{ sup E{t)) . 
te[o,T] te[o,T] te[o,T] 

FinaUy, this estimate together with equation (j6.24|) shows that 

from which it foUows that 

sup (hWll4 + IIPoJ"^WIl4) <Mo+(5 sup E{t) + CTP{ sup E{t)) . (6.26) 
te[o,T] te[0,T] t6[0,T] 

6.6. Estimates for curl^ v. As a result of the inequality (|6.26p and the identity curlj, v = 
curluoJ^^, by readjusting the constants if necessary we have that 

sup {\\cui%v{t)\\l + \\pod^cui%v{t)\\l) < Mo + S sup E{t) + CT P{ snp E{t)) . 
te[o,T] te[o,T] tG[0,T] 

Corollary 6.7. 

sup {\\ cmlr, v{t)\\l + WpoB'^ cml,^ v{t)\\l) < AIo + S sup E{t) + CTP{ sup E{t)) . 
tG[0,T] te[o,T] te[o,T] 

Proof. Letting act on the identity (|6.2p for curl^ w, we see that the highest-order term 
scales like 



D^cmluo+ [ D'^vDvAAdt'. 
Jo 



We integrate by parts to see that the highest-order contribution to curl^ v{t) can be 
written as 

curl uo- [ DS [Dv A A]tdt' + DS{t) Dv{t) A{t) A{t) , 
Jo 

which, according to ()6.26p . has L^(ri)-norm bounded by 

Mq{6) + 6 sup E{t) + CTP{ sup E{t)) , 
te[o,T] te[o,T] 

after readjusting the constants; thus, the inequality for the 77'^(ri)-norm of curl^u(t) is 
proved 

The same type of analysis works for the weighted estimate. After integration by parts in 
time, the highest-order term in the expression for poc^* curl^ w(<) scales like 

poB^ curl Uo - [ pad^Dr] [Dv A A]tdt' + poB'^Drjit) Dv{t) A{t) A{t) . 



Hence, the inequality (|6.8p shows that the weighted estimate holds as well. □ 

6.7. The a priori bound. Summing the inequalities provided by our energy estimates, 
the additional elliptic estimates, and the estimates for curl^ v shows that 

sup E(t) < Mq + CT Pi sup E{t)). 

tG[0,T] te[0,T] 

According to our polynomial-type inequality given in Section[431 by taking T > sufficiently 
small, we have the a priori bound 

sup Eit) < 2Mq . 
te[o,T] 
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